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Abstract 

The height probabilities for the recurrent configurations in the Abelian 
Sandpile Model on the square lattice have analytic expressions, in terms of 
multidimensional quadratures. At first, these quantities have been evaluated 
numerically with high accuracy, and conjectured to be certain cubic rational- 
coefficient polynomials in 7r _1 . Later their values have been determined by 
different methods. 

We revert to the direct derivation of these probabilities, by computing 
analytically the corresponding integrals. Yet another time, we confirm the 
predictions on the probabilities, and thus, as a corollary, the conjecture on 
the average height, (p) = 17/8. 
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1 Introduction 



The Abelian Sandpile Model is a non-equilibrium system, driven at a slow steady 
rate, with local threshold relaxation rules, which in the steady state shows re- 
laxation events, called avalanches, in bursts of a wide range of sizes and critical 
spatio-temporal correlations, obtained without fine-tuning of any control parame- 
ters. We refer to the introductory reviews [IH1]- 

In the set of stable configurations in the Abelian Sandpile Model on (portions 
of) a square lattice, at each site i e Z 2 , the height variable can take the values z% = 
0,1,2,30. Particles are added randomly and the addition of a particle increases 
the height at that site by one. If this height exceeds the critical value z c = 3, then 
the site topples. On a toppling event, its height decreases by 4 and the heights at 
each of its nearest neighbors increases by 1. 

A very natural question is: what is the asymptotic (i.e., infinite- volume) prob- 
ability Pi for the heights Zj, for i — 0, 1, 2, 3, in the ensemble of recurrent configu- 
rations? 

After some numerical studies [5HZ], the first exact result [8] concerns with the 
probability Pq for a site to be empty: 
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An analytic expression for the other probabilities was obtained in P|[TU]: 

1 3 2 12 It 

p ^2-T.-Ti + Vi + f< < lb > 



13 1 12 h 31 



4 2tt 7i z it 6 2 32 

Pa = I-^ + ^ + T + 32 ; < ld > 

where I\ and I2 are expressed as multiple integrals. These results had been ob- 
tained by using a mapping from the set of recurrent configurations onto the set 
of spanning trees covering the lattice. These trees are rooted (on the boundary 
of the lattice, where dissipation occurs). Introduce the concept of predecessor: a 
vertex j precedes a vertex i if the unique path on the spanning tree from j to the 
root includes i. Then, the probabilities Pk at the vertex i are simply related to 
the numbers Xk of spanning trees in which the vertex i has exactly k predecessors 
among its nearest neighbours. And the X^s can be, at the end, expressed in terms 
of the lattice Green function. 



1 Some authors prefer the values z\ = 1, 2, 3, 4. Results are easily translated between the two 
notations. 
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Furthermore, an indirect argument fixes a relation between I\ and I2 (see equa- 
tion (I15p later on)H Thus any single further linearly- independent information on 
l\ and I2, or on the P^'s, would have fixed the height probabilities completely. 

An extensive account on the derivation of these results is provided in [H], 
together with several other interesting properties^. 

As a corollary, the average density in the ensemble of recurrent configurations 
is given by 

/ * x — v -. ^ 7 3 3 37o . . 

(p )=J2kP k = - + - j + 2 

' 4 2tt tt z 32 

k=0 

As reported in [TJ, this quantity was conjectured by Grassberger to be@ 

(P) = ^ ■ (3) 

An interesting observation is the following: the expectation value of the height 
probabilities don't help in the understanding of the conformal features of the cor- 
responding field theory in the continuum, at least in the whole plane. But it 
is not so in the presence of a boundary. Indeed, in [HHI1], the evaluation of 
the height probabilities in the upper half plane has been used to reveal that the 
continuum theory is a logarithmic conformal theory with central charge c = —2. 
Afterwards, also two-point correlation functions for the height variables have been 
computed [15] and found in agreement with the prediction of a logarithmic con- 
formal field theory based on field identifications obtained previously. 

A rejuvination of the interest in the exact determination of (p) has arised 
with the work of Fey, Levine and Wilson [15] [T7], in which a subtle difference 
has been elucidated between the uniform average on the ensemble of recurrent 
configurations, and the properties of the critical system with conserved mass. As 
the discrepances in densities between the two regimes are numerically very small 
(Sp/p ~ 10~ 4 ), although the numerical determination of the integral I2 appearing 
in (J2J) has a much higher precision, it would have been more satisfactory to have 
an exact result for at least one of these two quantities. 

In fact, it was possibly in part this rejuvinated interest that led some time 
later to two independent proofs, methodologically similar, of the density conjecture 
(and, through the argument above, of all the height probabilities) [18][T9]. The 

2 A certain combinatorial quantity, know to be finite, is formulated as a lattice integral present- 
ing a divergence: tuning to zero an overall factor in the divergent part gives the forementioned 
relation. 

3 In [14], in equation (4.1), the authors also correct a misprint in equation (32) of [10] which 
was wrong by a factor 2. Here we notice a misprint in their equation (4.8), where the second 
term on the right-hand side should be 7/w instead of 7/ir 2 . 

4 Thcsc authors use the range 1 < Zi < 4, and accordingly write (p) — 25/8. 
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single missing linear relation has been the intensity of the loop-erased random walk 
at a first neighbour of the source of the walk, that is combinatorially related to 
the density (and turned out to be 5/16 on the square lattice). The role of the 
loop-erased random walk in the uniform spanning tree model (and thus in the 
Abelian Sandpile Model) should not be surprising, since the works on the subject 
culminating in Propp and Wilson exact sampling algorithm [20J. 

In this paper we shall provide a different proof, conceptually simpler (although, 
admittedly, theoretically less illuminating): we shall revert to the original formu- 
lation of the problem, and evaluate exactly the integrals in question. They have 
the form of two-loop Feynman integrals in a two-dimensional scalar field theory 
on the lattice. There has been a long-time effort in order to reduce the 

evaluation of lattice Feynman integrals at one-loop level, through simple al- 
gebraic methods, both in momentum space [2TJ|22], an< ^ * n coordinate space [23] . 
These methods have found also important applications in two dimensions, respec- 
tively [24"tl2o'] and [25] . In particular, in [27] there is an extension to the triangular 
lattice, which can be of help to generalize our procedure to this system (remark 
that the study of sandpiles on the triangular and honeycomb lattices, with the 
aim of height probabilities and correlation functions, has also been considered 
in [321E8]). 



1.1 The integrals to evaluate 

We shall use the following notations of lattice momenta, which are common in 
lattice field theory 

p^:=sinp M (4) 
P» ■= 2 sin -± (5) 

where, in our two-dimensional case, the index \x can take two values that we choose 
to be and 1. Then 

f ■■= E pi ( 6 ) 

m=o,i 

is the quantity, invariant under the lattice symmetry, which appears in the lattice 
propagator 

p z + h 

We have added the regulator of the infrared singularity h just to have well-defined 
quantities in all our manipulations, but we are interested only in the limit of 
vanishing h (and thus, to integrals for which this limit exists). In the following we 
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will not mention explicitly the regulator h, and the extraction of the limit will be 
understood where pertinent. Given the shorthands 



/ dp := £ tr £ t 1 ; 62{p) := (27r)25 °° o) 6{pi) ; (8) 

we want to integrate polynomial expressions in the lattice momenta (Jl]), in the 
measure 

d/i := dpdqdk5 2 (p + q + k) A(p)A(q)A(k) (9) 

which is invariant under all the permutations of the momenta p, q, k, under ex- 
change of the indices with 1, and under simultaneous inversion of all the mo- 
menta along one of the lattice axis. These invariances imply relations between the 
integral of different polynomials, to which we will refer generically in the follow- 
ing as "symmetry of the integration measure". In particular, we use the symbol 
A — > B to denote the fact j d[iA = J dfiB. 

In order to define the integrands pertinent to the expressions in ([1]), we have 
to start form the matrix M(ci, 02,03), given in [T4J eq. (3.18)] 



Cs e «po+«go e -iqo+m e -«po 
Ca e ipi+iqi 1 e~ ipi 

-ipi-iqx p 2iq\ ipx 



(10) 



/ci 1 e iqi 1 \ 

M(c 1 ,c 2 ,c 3 ) = 

\c 2 e 

The interesting quantity is the integral [HI eq. (3.17)] 

J(ci,e 2 ,c 3 ) = J dfj,iship det M(c 1 ,c 2 ,c 3 ) . (11) 

It is soon realized that the integral is real, does not depend from C3, and is, of 
course, linear in c\ and c 2 : 

I(ci,c 2 ,c 3 ) = \ (J\Ci + J 2 c 2 ) . (12) 

o 

The factor 1/8 is due to our choice to maintain the usual definition of the lattice 
propagator. This differs from the choice in [9"|llUpi4j by a factor 2, and we have three 
propagators in the definition of the integration measure. Then, for the quantities 
Ji, I 2 defined above in ([T|), 

h = Ji+^-iJ J 2 ; (13) 

16 4 

h = 8 J 2 - — + — . (14) 
7r ir z 
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In [H] it is shown, by an indirect compatibility argument, that the relation 



2 4 

Ji + h = 5 (15) 

must hold. This is verified numerically to high precision (to order 10~ 12 ). It is 
also observed numerically to high precision that 

J2 = \ ■ (16) 

If these relations hold exactly then 

1 1 3 12 , . 

Pl = Z " 9 2 + Z5 ; ( 17a ) 

A = - + 3 ; (17b) 

8 7T 7T ,} 



3 114 . . 

p 3=o-^ + ^ + ^; (17c) 



2n TT 2 7T 3 

and also the conjecture by Grassberger on the density follows 

fc=0 



1.2 The strategy 

The details of the derivation are given in the following sections. Let us outline 
here the general strategy that we adopted all along the calculation. By exploiting 
the symmetry of the integration measure, we try to obtain, in the numerator, a 
factor which can cancel one of the propagators (at vanishing regularisation). Say 
that we get a p 2 in the numerator. Then, we write all other appearences of p^s as 
— {% + &fi)'s. Thus, possibly through a trigonometric expansion at the numerator, 
the remaining integrals are factorized in independent one-loop integrals in the two 
other momenta. Some useful trigonometric identities used at this aim are 

ju=0,l M=0,1 

K = €+pI~\€pI + 2 ^- ( 2 °) 

The latter, which is valid when the sum of the momenta p, q and k vanishes, is 
sometimes useful also in the inverse form, in which q^p^ is expressed in terms of 
the rest. 
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We shall need the very elementary one-loop integrals: 

<« = 2; /*p = o ; /*f = 5i = < 21 > 

From these building blocks, other integrals soon follow, for example 

# = /dg ^ 2 -^ -!; (22) 
-2 /• / ^2 



dq T = J dq \¥ + 4¥) = 2-^ ] (23) 



and the slightly more tricky 



dq = - J dq — = - / dq q q x = 2 . (24) 



We shall also need 



go 2 cosgi [ (ql \ql lggg 1 q%$\ 1 1. 

g 2 7 ^ V? 2 4<f 2 <f 8 q 2 J ' 4 tt 



One more trigonometric identity has been used in the Appendix in order to com- 
pute a slightly more complex integral. 

2 The integral J 2 

The contribution proportional to C2 comes from the integral of 

2 sin Pq sin qi {sin p 1 [( — cos g + cos <?i ) sin p — ( 1 + cos p + cos Qi ) sin <Zo] 

+ singi[(l + cospi + cos go) sinp + (cospo — cospi) sing ]} (26) 

which has been obtained, from the definitions (TT0l - [T2|) . by performing a sym- 
metrization by first changing all the signs of Po,qo, and afterwards of P\,q\- A 
convenient rewriting, using k + p + q = 0, is 

2pb"gT [(poW- %Pi) + (%h - k -qj) + (& Pi -pb"&i)] , (27) 

that is, through a complete symmetrisation, 

X 2 

2 [{PoQi- QoPi) + (Qoh - k ~qj) + (fcopT-po k ^)] ■ ( 28 ) 
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The expression in square brackets has a nice geometrical interpretation: it corre- 
sponds to twice the area of the triangle with vertices located at the points (p~o, Pi), 
(qS, gl) and (k , fa). 

A different rewriting, reverting to p and q only, and changing p with q in some 
terms of (1271) . is 



(P0?i-?oPi) + 2poPo + <?o<?i (qi-Pi) + 2piPi + qiqo (ffo-Po) • (29) 

The last two summands clearly give identical result after integration, and we have 
to evaluate two integrals 



J. 



(a) 



J, 



so that J 2 = 4 a) + 4 b) . 



(b) 



(pogi - go pi) 2 ; 

dfj,Ap^p + q -qj (qj-pl) 



2.1 The integral J 2 (a) 

We start by evaluating the integral in f[3"Uj) . We rewrite it as 



(30) 
(31) 



• M=0,1 ' ^=0,1 ' V=0,1 

In the integration, exploiting the symmetries, 



/i=0,l 



Is 2 



iy=0,l 



P 



4 J m 

M=0,1 



s2-2 -2-4 



i/=0,l 



/i=0,l 7 v i/=0,l 



(32) 



(33) 



Using (|20|) in the inverse form, the expansion of the square produces 3 contribu- 
tions. The first one is 



-2 -2 
, p 



^ M=0,1 M=0,1 7 



(34) 



-2-2 , d 7*2-2-2 

-TP + TE k 9 P 
lb 
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where we used also the fact that in the integration the factor k 2 cancels all the 
dependence from k and the subsequent integration of vanishes, while the 
integration of q 2 u p 2 does not depend on the values of fi and v. 



The second contribution is 



hk 2 -q 2 -p 2 ) qlPl^-\k 2 q 2 p 2 + \p 2 %(q + k) 



-j> 



1 f2«2«2 -2 -4 i 1 f2-2-2 ^ -2 7*2-4 

- -k q p + p q + -k q p - -p k q . 



We now combine the last term with the similar one in (1331) . that is 

2n i 



1 

16 



i2"\ -21 2 



16 



2,;2 _^2-2l2 



.4=0,1 ' ^zv=0,l ' ^M=0,1 

ipg 2 (Pg 2 -^ 2 ) 

m=o,i 

1 7*2-2-4 1 7*2-4 -2-2 ^ 7* 2 -2 -2 -2 , ^ 7*2-2-4 -2 



— » 



M=0,1 



16 



16 



(35) 



(36) 



. 1 ; 2 -2 -4 1 f 2 -2 -2 -2 

-> p q - — k p q q . 
By collecting all the pieces, and using the elementary integrals fT2"Tj) . we get 

*w + (I - 1) * w ■ (37) 



In conclusion, we find that the value of the first term is 



J 2 ] = f dfi (p q 1 - g Pi) 2 = ^ - . 



(3? 



2.2 The integral J. 



(6) 



We now consider the evaluation of the integral in (J3TJ), and remark that the inte- 
grand can be written as 



4p Po + go qi (?i - Pi) -> -2 ^Ppkp V" ~~ ?0 

^=0,1 i/=0,l 



(39) 



because the terms with /i = v give a vanishing contribution (they are anti- 
symmetric under the exchange of p with q). Repeated use of ( TT9|) and ( 120]) gives 



n=o,i 



u=0,l 



u=0A 



■ (40) 



We split this evaluation into three terms (the following (jUJ), (T42]) and f )45|) ). A 
first contribution is 



/i=0,l 

{q 2 -k 2 -f 



Ck 2 -3q 2 -p 2 



2 ^2 pIK 

li=0,l 



q 2 = -3q 2 J2p,K^ 



(41) 



Ai=0,l 



which vanishes in the integral. A second contribution is 



q 2 — k — p 2 
4 



u=0,l u=0,l 



i/=0,l 



^=0,1 



We have now, in all summands, an exposed propagator. Following our general 
strategy, we rewrite the remaining expressions using p + q + k = 0, namely 



-J 2 £^ = k £ [#(i-i#)+#(i-J#)+2* 



!/=0,l 



!/ = 0,l 



*vPv 



-+ U 2 £ fef 1 -^) +3^(i-i^)e(i-i* 



i/=Q,l 



(43) 



— >■ 



ig 2 £ (% + 3klpl-2ktpl + fk 



4 -4 



i/=0.1 



so that the whole contribution from (1421) is 



« 2 £ 

i/=Q,l 



| - ktPl + I ktPt) -+ I q 2 k 2 p 2 - q 2 p 2 kt + \ q%f 



We are left with the third term 

I J2pIK £#(# + 



0=0,1 



u=0,l 



One summand gives 

I E pIK £ ^ z^ 2 (% 4 + ^) 



(44) 



(45) 



(46) 



m=o,i 



v=Q,l 
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(47) 



because 

dfitfo k lq\ = I d ^Pl k l (<f - & = J d ^Pl k l qt - Zqlq 2 + (q 2 ) 2 

dfxpl k 2 + q 2 (qf - qf)] = J d[ip\ k 2 Q q A Q 
where we used the fact that 

d/iplk 2 Q q 2 (ql -ql) = [ dpdkp 2 k 2 fop\ - \p\k\ - 2p 2 + i#^A (p)A(k) 



4 1/ 4 
2 2 

7T 2 V 7T 



4 8 



7T 7T Z 



0. 



The second summand is 
1 

8 



! £.2 \^ -,2 .2 , ' . 1 12 -2 , 1 ^2 ,.2 ,2 .-,2 



M=0,1 



i/=0,l 



Pv -> + 7^0 fc 0Pl ?1 



(48) 
(49) 



and 

while 

and 

so that 



^0 ftgp? <? 2 -> ^0 k lPl (<f - ql) = 1(2 - \ 



1 ^2 -2 T 2 -2 . 1 ~4 r*2 ~2 1 -2 -2 F 2 -2 

^p Pik qo -> ^PoMo - ip pqK% 



plPlKql (50) 



(51) 



- / d^p P A; g = -- / dfcdg 



A- 2 « 2 
I 3 2 



~ + IS (52) 



7T 7T Z 



/x=0,l k=0,1 

By collecting all the pieces 



3 2 
4 2^ ^2 



(53) 



(54) 



- fdnJ^plK ^2 €(€+$,)= I d ^Pl k lql + \ 

■> m=0,1 1^=0,1 
and using the result (1751) . computed in the appendix, the whole expression ^ is 

-2/4. 

M=0,1 i/=0,l 

^ (*J *g ql + I q 2 kY - q 2 P% + \ q%pi) + ~ - 2! + ^ ( 55 ) 

6 6\54413 2 11 
+ Z ~ ^2 J + g ~ n + ^ + 4 ~ 2^ + ^2 ~ ~8 + 2^ ' 
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In conclusion, by adding the first and the second result, computed respectively in 
(J38D and (J55D, we get 

^ = A_±_I + ± = ± (56) 

8 16 2vr 8 2tt 16 v ; 

in agreement with the prediction J 2 = 1/2. 



3 The integral J\ 

As anticipated in the introduction (see equation f[T^]) ). it is expected (by an indirect 
argument) that 

- - - — (57) 

Similarly to our evaluation of J 2 (but, as we will see, in a simpler way), we can 
attack directly the evaluation of J\ + J 2 , and produce an independent check of 
the relation above. Recall that J\ is the contribution to (TTTj) proportional to c\, 
namely 

Y = J d/i2p^'qj(ko qi - Pi + gb>i - h + pb~k 1 - qi) . (58) 
Writing p — q = pcosq — qcosp, restate the integrand above as 

2pb~gT [(%Pi-PoQi) cosfci + (pSk 1 - kopt) cosgi + {k ~qj - h'qt) cospi] (59) 
and remark that, by the replacing 

cos9 = l-^9 2 (60) 
all the contributions in which we take the 1, that is 

2po qi [(3o Pi - Po qi) + (Po fa - k pi) + (k qi~fa go")] (61) 
are exactly — J 2 /8, thus, if we keep only the other terms, we have 

j _|_ j r p A 

8 2 = ~ / dfip^qi (^pl-p^W)kl + (p^fa-k pl)q 2 + (koqi-fa^pl . 

(62) 

Manipulate the integrand by exchanging q with p, and index with 1, to get 



-Po qi 
In conclusion 



(go Pi -po?i) \k 2 + (po fa - hp^f 



-> Po 2 q\ \k 2 + po 2 cos p\ ki q 2 . 

(63) 



1 , , , x 1/1 1\ 2 /l 1\ / 1 1\ 1/1 1 



as it was predicted. 
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4 Conclusion 



We have been able to analytically compute some lattice integrals that, through 
the work of [HUffi! EH E] and references therein, describe the height probabilities 
in the ensemble of recurrent configurations of the Abelian Sandpile Model on the 
square lattice, in the thermodynamic limit. 

The numerical values of these integrals were already known with high precision, 
and the exact expressions solidly conjectured, as rational-coefficient polynomials 
in 7T _1 . Most importantly, a recent indirect calculation of statistical properties of 
the loop erased random walk, or equivalently of domino tilings with prescribed 
local patterns of monomers and dimers, was sufficient to determine completely 
these values [T8UT9]. 

Nonetheless, our direct evaluation of the original lattice integrals, with their 
strikingly simple results, could be of some interest, and of some use for future work 
in similar contexts. 

Let us stress again that this result is not based on any new deeper understand- 
ing of the properties of the sandpile model, but completely relies on elementary 
trigonometry and symmetry considerations, mainly with the aim of reducing two- 
loop lattice integrals to quadratic polynomials in one-loop integrals. In particular, 
at some point we used results previously obtained in [21]. In principle, we do 
not see any obstacle to recover similar results on other two-dimensional regular 
lattices. 

A One more integral 

We need the evaluation of the integral 




(65) 



We first observe that 



(po go) h 




(66) 



but it is also 



(po ko) (go ko) = -(q 2 - p 2 Q - fcjj) + -p 2 k 2 -(pjj - q 2 - fcj}) + -ggfcjj 



(67) 




(68) 
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so that 

HpUl 6# - 12 PUl + Aplqlkl + 2p*q 2 - f . (69) 
Let us start with 



6 J dfi (p$ - 2p 2 q 2 ) = -12 J dfj, (2^ 2 + ApSq^) = -12 G 2 = ~ 



(70) 



(72) 



where the integral G 2 was defined and calculated numerically in [29] and subse- 
quently computed in [2J1 eq. (A. 9)]. Then, 

4 f d^g$g*g = 4A< 8 > = ± ; (71) 

where the integral had been introduced and computed in [2U eq. (A. 6)]. For 
the evaluation of the last term we use the same trick that was used in [21] to 
compute A^ 3 ', that is use the fact that 

J dfip 4 (2q 2 -pl) = Jdfi (2piql-pt) 

= jdfi(p 2 -pl) 2 [2q 2 -p 2 -(2q 2 -p 2 Q 
which has the consequence that 

J d»pt (2 q 2 -Pl)=ljd^ [P\{2 q 2 - P 2 ) - P 2 {P\ - Pl){2 q 2 - p^)] 
= \ Jd^p 2 [2qi-pt- (Pj - Pl){2 q 2 - p 2 )] ; 
so that we are left only with elementary evaluations, that bring us to 

fd^ {2ql-P 2 )=-\ + --- 2 . (74) 

J 4 71 7T Z 



(73) 



In conclusion 



,::.,.> 1 1 5 6 6 ,6 6 
dfikfaftfi = -- + --- + = -1 + . 75 

4 7T 7T 7T 7T Z 
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